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ABSTRACT 


A  basic  conical  element  subjected  to  axisymmetric  linearly  varying 
teiqperatures  on  the  inner  and  outer  surfaces  is  solved  for  internal 
stress  and  compatible  ed^e  loading  and  defoimation.  Uniqueness  of 
the  result  is  proven. 


Descriptors 

Cones 

Thexmal  Stress 

Elasticity 

Shells 


1. 


N0TATIC9I 


T 

AT 

T 

s 

X 


✓ur 


KJ 


s  ) 


oc 

V 

E 

b 


K- 


Et" 


Temperature 

Temperature  difference 

Mean  temperature 

Coordinate  distance  from  apex 

Half  opening  angle  of  cone 

TOiermal  constants 

Meridian  angle  change 

Normal  displacement 

Axial  stress  restiltants  in  neridional 

and  circumferential  directions 

Stress  moment  resultants  in 

meridional  and  circumferential 

directions 

Coefficient  of  thermal  expansion 
Poisson's  ratio 
Young's  modulus 
Thickness  of  cone 


(degrees  Eahrenheit) 
(  "  "  ) 

(  "  "  ) 

(in.) 

(degrees) 

(degrees  Fahrenheit) 
(dimensionless) 

(in.) 

(Ih/in.) 


(Dh.in./ln) 
(ln/in/°F) 
(dimensionless ) 
(Th./in"^) 

(in.) 


Bending  rigidity 


(lb. in.) 


-  Principal  xadius  of  curvature  in 
circumferential  direction 
c|)  •  Angle  between  shell  normal  and 

cone  axis 


(in.) 

(degrees) 


2. 


KOTATICai 


A.  s  Qj  -7— 

(j  _  - __i — _  Stress  function 


(Ib./in*) 


(,,.y= 


A. 

A-s 


(... ) 


DC  ss  '3  j  >s 
Etr  . 


>“  = 


>  Coordinate  transformation 
-  Parameter 

>  Constants 

L(...)=  SC.,.)"  +  (..,y  -  r  (...)  s  VLc...) 


L  (  )=  +  i-  (..  ,V  -  1  c...) 

C: 

X 


X  X 

Constant  of  integration  c  -  •  j  •  •  >  B 

Horizontal  shear 


(Ih  ./in.'^) 


SOLUTION  SUMMABI 


r 


r  AT 

=  T„-T; 

[  T 

_  T„  +  T; 

X 

(  AS 

=  s^-s, 

A  - 


c  = 


X  -  - 


Thermal  Cknsstaats 

t^-t' 

As 

t  A  s 

Edge  Loads 

Oi  VC  ( )  c 


Edge  deformations 


X.  -  —  cx  A  s  ^ 


r 


6  = 


D  - 


s.t"' 

As 


5.AT'-s,  AT^ 


t  As 


^  0(  \C,  C  \  V'  )  \^  A  fevu  ^  'V'  C  S  "V  U 


vXO-' 


^  -^...A.  sXt;. 

A 


A 


Internal  Stress  iresultants 
—  •VSJ^  (XKCI'V'V/)  G  ^ 

_  a  l<  c A  +  Cs -v-  D  ) 

*>  t/  1 


it. 


IHTRCDOCTION 


[  The  solution  of  thermal  stress  in  an  aihltrary  thin  shell  of  revolutian 
subjected  to  an  axlsytametric  temperature  distribution  specified  on  the 
inner  and  outer  surfaces  of  the  shell  may  be  solved  by  the  cone  coupling 
technique  developed  in  SM-385OOJ  "Analysis  of  a  Shell  by  the  Truncated 
Cone  Approximation"  by  M.  B.  Harmon,  provided  that  a  rigorous  solution 
of  an  individual  conical  eleinent  is  known. |  It  is  the  purpose  of  this 
report  to  develop  such  a  basic  solution  for  a  linear  ten^jerature  distri¬ 
bution  defined  at  the  edges  of  a  truncated  conical  element. 

The  governing  equations  are  the  well-known  Beissner-Meissner  shell 
equations  as  extended  by  Melan  and  Farkus  for  thermal  loading.  These 
consist  of  a  coupled  set  of  second  order  ordinary  linear  differential 
equations.  Particular  integrals  are  developed  by  means  of  the  method 
of  undetermined  coefficients  in  a  polynomieJ.  for  each  dependent  variable. 
It  is  shown  that  the  particular  Integrals  are  unique  and  represent  the 
desired  basic  solution  provided  that  compatible  edge  shears  and  monents 
are  provided  at  the  edges  of  the  conical  frustium. 


5. 


CaEOMETBY 


The  notation  employed  will  be  that  of  Fluegge  to  be  consletent  with 
previous  developments  by  the  author  [3]  ^  [S']  . 


Co>1^  VHici  iOlT'f'CvCe. 


roffEBATUIjE  DESCEIFTIOW 


In  1611118  of  the  nonoal  coordinate  z  to  the  shell  mid-surface,  measured 
positive  outward,  a  Taylor  expansion  of  the  temperature,  T,  ebave  scaae 
arbitrary  stress-free  datum  would  be, 

T  =  r^cs)  X  r.  u)  +  ^  ^  < 

^  ^  ~  thickness. 

Now  since  "t «  S  ,  one  may  ignore  the  higher  order  terms  and  write : 

T  =  T.csl  +xTu)  (1) 

This  linearization  is  consistent  with  the  approximations  inherent  in  the 
Love-Kirchhoff  formulation  of  thin  shell  theory. 


SHELL  EQUATIQMS 


The  Reissner-Melssner  equations  in  terms  of  T^  and  T^  are  [ll  ; 


-+-  S  X  -  X  =  ocC'-^^)  S^T  —  (2) 

1  K 

s'^U-^sU-U  =  (sc^p  X+ oes^T^ ')  -  (3) 


Where: 


X 

K 


cLlaj' 

d'S 

Ix  Ci-v^) 


Meridian  angle  change 


=  Bending  rigidity 


o4  =  Coefficient  of  thermal  e^qpanBion 


H.  =  Young’s  modulus 


=  Poisson's  ratio 


u 


-  ^  6  Q: 


t’ 


^ ®  SiiJTSSS  ^Vuxci/ jLo£t 


jL 

^Ls 


Note  that  this  definition  of  U  differs  from  that  of  FLuegge  by  a  factor 

of  .  This  unessential  modification  is  due  to  Melan  and  Baxtais  and 
t^ 

involves  the  particular  integral  only. 


8. 


lEMrEEATUEE  CONSlEAiraS 


Let  the  specified  temperature  dependence  upon  s  be  given  by  1die  following 
linear  e3q)reBSions : 


=  A  5  B 
Tj  -  C  s  +  D 


ih) 


The  unknown  constants,  A  B  C  D,  may  be  evaluated  from  temperatures 
defined  at  the  outer  and  inner  surfaces  at  points  s  =  s^^  and  s  =  Sg. 


t 


5  —  "S  ( 


T a  -  As,-fB  +  '^  (Cs,4-0) 


r':  ^  As.  +  B  -  4  CCs.i-  0) 


5  =  s- 


"  A  s-^.+ B  ^  (C  3^+ d) 


^  T”^'  :=As^+B-  —  (Cs^+O) 


outer  tenperature 


inner  ten5)erature 


outer  temperature 


inner  taoperature 


9- 


These  are  four  e<iuatlons  in  four  unknowns. 


1 

5.  ! 

A 

1—. 

i 

1  o 

-A  -1 

6 

\ 

T: 

S  1 

1 

— 

•  0 

1 

-1 

i  D 

T-: 

C 

with  solutions  given  by: 


A 


A  s 


-  s.T 
As 


at"'-  at' 

t  A  5 


D 


•t  A  5 


where 

AT  ^  r,-T; 

=  Ten5)eratuie  difference 

=f  _  T+r; 

2 

=  Average  teii5)erature 

As  =  s^-  s, 

=  Slant  height 

10. 


GOVERHING  EQUATKBS 


Differentiating  equations  (4), 

t  =  A  ,  f ,  =  C 

and  substituting  these  values  into  the  Reissner-Meissner  equations 
(2),  (3)  one  obtains. 


s'"x!.  +  sx-)c  =  s  U  +  occ.i+v')  s^C 


4E 


5U+SU-U  =  — (  s  X  ^  K  ') 


For  convenience  of  manipulation,  define  the  constants. 


pr  _ 


H  - 


4^^ 


6,*(3 


Or  -  (l+'O  C 


A 


(6) 


■T-  _  ^  H  1 

X  -  - ; -  Oi  A 


The  governing  equations  become: 


-X  -  - 

~  V"  5  (J  -+•  Gr  5^  — 

—  (7) 

IJ  ^  s  \J 

-  0  ^ 

M  s  X  +  X  6'" 

(8) 

This  system  is  a  set  of  two  ordinary  coupled  second-order  linear  non- 
homoepjneous  differential  eqiuatlons  vhose  complete  primative  consists 
of  a  fundaioental  set  of  solutions  associated  with  the  hoosogeneous  part 
and  a  particular  integral  of  the  non-homogeneous  equations.  P* 


HCMOGENEOUS  SOIiJTICHJ 


The  homogenecns  equations  are: 


=-F5U 


s'"U-+'SL)-U  -  MsX 


Define  the  cone  operator  [2^  p.  196, 


LC..)  =  s(.^,V\(.,,V-4:  , 


then 


LCx)  -  -FU 


L(U  )  ■=  H  X 


(9) 

(10) 


(11) 

(12) 


12. 


Successive  opei*ations  yield, 


i.e. 


LL  (X)  =  -  FL(U) 
L  L  (o'!  =  Hi.  (Xl 

LL  (X)  +  FH  X 
LL  CU  )  +-  FW  U 


=  -  F  H  X 
=  - F  H  <0 


-  a 


Now  U 

be  written: 


4sGls 

V 


so  that  the  second  equation  may 


LLCsi^c,]  +  FHsQc  =  o 


W 


In  order  to  solve  this  equation,  one  may  factorize  the  operator  and 
obtain  the  two  equations: 


1 _  (_-SQ^')  i  L~X  "  o 


(15) 


where 


i<  r 


(16) 


13* 


The  solution  of  eqiuatlons  (15)  laay  te  ctotained  in  terms  of  the  Kelvin 
function  [2*^  ,  ber  x,  bei  x,  ker  x,  kei  x  and  their  deriTOtives 
employing  the  coordinate  taransformation, 


oc  =  Ts  -  (IT) 

*  The  Thompson  functions  ^oted  by  Iluegge  are  the  Kelvin  function 
since  Kiongison  was  also  known  as  Lord  Kelvin. 

to  give^ 

sQ,  =  +  qCbe.-x,  +i  b^r'x)  'I 

Vlhere  '  '<  0^  constants  of  integration . 


An  identical  solution  for  X.  may  be  written  down  in  terms  of  new 
constants  of  integration  ,  Cq‘,  however,  the  two  sets  of  constants 
are  related  by  equations  (ll)  and  (12)  and  it  is  easier  to  cibtaln  the 
solution  of  X,  directly  finm  these  relations. 

Define  the  fundamental  set  of  solutions  by  equation  (I8), 


j 

-  he  ir  X,  ■ 

_ 

■% 

lo  Cl.  ^70  ^ 

-  lo  e  t  -pc 

'L 

X 

heAx, 

^3 

—  K  e-  V"  tl.  ' 

.  1. 
X. 

~  \Ctl  X- 

X 

Kev-'x 

Then  the  homogeneous  solution  becomes  briefly. 


ih. 


Now  from  egufttion  (L2), 


X  = 


H 


Hu) 


t 


4  £  (Uyt*  “t 


LlsQJ 


X,  =  LCsQj 


s  = 


on  •  - 

•  ^4 

we  must 

define  L 

in  terms  of  x. 

zc 

■  Ji 

A 

dx. 

X  / 

A.  ) 

(-  ' 

A_ 

4-  -L- 

X 

Lx,  1 

cLx. 

X 

operator 

L  ( 

..  ) 

from, 

(..;)  = 

s  (.. 

(,. 

.)'- 

I  c.. 

,) 

4(' 

/ 

.X  - 

4 

c... 

)]  = 

-  > 

"lC.. 

) 


, ) 


where  (  .  .  .  )' 


± 

cLx, 


The  expreesion  for  X  becomes: 


15. 


Due  to  the  linear  property  of  the  operator,  one  vrites: 


L  ) 

=  L,  (.  bev'x  -  — 

be:^x)  := 

-loeoc--  be^'x  := 

L(4xJ 

=  L  ^  be:  X  +  ^ 

bet^'x)  = 

bev'^.  be: So  =  (20)  ^ 

L 

LC^^e^x-" 

ket  X  )  - 

'  be:  x  -  -  ke^'x  -  =•  - 

L  (^4) 

-  L  0<e:  X  -h  - 

X 

kew'x)  - 

b  e.  X  -  -^  VCdl  X  -  s  Au  j 

X 


fet 


io 


(21) 


Where  one  has  used  the  relations  [2]  p.  171 


-7.  -  -  be 


'  L  ' 

~  X. 

J 


b  e  u  X 


■=  - -  bcu-jc 


V.t.v  X- 


7w 


>y  ’ 

vec  L  7U 


Ke^x  l^cc'x. 

X 


16. 


PARTICULAR  INTEGRAL 


Turning  now  to  the  particular  integral  of  (7)  (8)  and  guided  by  the 
polynoaieQ.  form  of  the  coefficients;  a  solution  by  the  method  of  un¬ 
determined  coefficients  is  tried. 


X-  -  Cq  -t  C,  S  -V 

U  =  Of,  +  0,  s  -v-  Dx 


(22) 


c,  C.  and  Do  D,  are  not  related  to  previously  defined 

constants . 

Substituting  equations  (22)  into  equation  (t)>  one  obtains, 

-h  £  -(Cf,  +  6.s  Fs  (D.+D.i  +  O^S^') 

-Crs^  =  O 


Equating  coefficients  of  s. 


(o)  -  ^  o 

(5)  C^-C,  t-FDf,  =  O 

{s^]  0-  +-  - 

(s")  F  D,_  -  o 


C„  ^  0 

D  c)  -  O 

3  F  0,  -  C-  =  o 


O 


17. 


Substituting  equations  (22)  into  equation  (8) 


+  £  CC'^+-aDi.s)  -  I.D,  +  D,s+0^s^')-^s(,C^  +  C^£-vC^s^) 


Equating  coefficients  of  x, 


-  T.  s'"  ~  o 


0=0 


Do  -  O 


0,  -D,  -HC,  =o 


C.,  =  O 


a  D,  ^10,  -  D,  -wc^-l  -  o  3  0^-HC,-T  =  o 


HC^  ^  o 


C^=  o 


The  final  solution  is: 


=  C,  -  D.,  D.  ■=  O 


4  (>+  V  )  C 

is'^toT  A 


=  -  oL  A  iXv-  (G 


X  =c^s  -  -c^As  ^ 


0= 


D  s  = 


4  (,  \  t- v)  c 

t"-  Ci>^  (6 


Q.  - 


18. 


One  may  now  write  the  cong^lete  primatives  of  (t)  and  (8) 


X  =  + +  -(25) 


vSQ^  '^4  -V-0L.1CS  C*"'"''^')  C-  (26) 


C<»islder  now  a  cone  closed  at  the  apex  and  extending  to  infinity  at 
the  open  end. 


When  the  shell  is  closed  at  the  apex,  one  obtains  the  boundary  conditions 
o  j  [23  p.  198#  while  and  Cg  are  associated  with  stresses 
which  are  non-zero  only  In  a  neij^orhood  of  the  bottom  edge,  i.e.  in  a 
neighborhood  of  Infinity. 


One  concludes  that  for  values  of  s  defined  in  the  open  interval 

O  <  3  <  =<> 

The  solutions  are: 

9C-  — '  *  S 

3  ~  J.  \<  6  C.^ 


19- 


UHIQUEIKESS 


It  follows,  moreover,  that  these  particular  integrals  are  unique. 

The  proof  is  comparatively  simple. 

I 

I.  The  cooiplete  primatives  aro  unique  since  the  coefficients  in  the 
differential  eijuatlons  (7)  and  (8)  are  analytic  [6]  p.h8  and  69. 

II . 

The  form  of  two  primitives  with  different  particular  Integrskls  4, 
4>^  would  be : 


UQs)^  =  f 


Where  ^  and  ^  are  the  homogeneous  solutions. 


Now  since 


^  o 


o 


'1. 


(sQ),  s  , 


while 


in  the  interval  0  <  S  <  oC 


one  obtains  necessarily, 


4>,  =  4>. 


20. 


STRESSES  ARD  DISPLACEMEHTS 


From  ex')  one  may  now  con^nite  the  remaining  forces  and  displacements. 

C  —  (28) 

--oCicCK^')  C  —(29) 

[x  -\'  j  >c  - U+v)  T,  ”1 

—  ^"U  +  '^)  o^A  /fclw  |2>  -  0)  3 

( 

^oL\<04j)|_At;uvja4-Cs  +  D^  —  (30) 

--t  I  T  -oJti+v)r,  ] 

=  [_-^\+^)^k  TZ^(2  -  o(U+^)  (Cs  +  D)] 

I^-v-Cs  +  D”\  —  —  (31) 

s 

AaT  -  -  {^  'X-  ct  ^ 

o 


/A  Aa-'w  ^  s 


^  .s'" 


(32) 
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CQKICAL  mJSTRUM 


If  a  conical  frustrum  "free  body"  is  isolated  from  the  Infinitely  long 
cone  and  loaded  by  reactive  edge  loads  of  magnitude, 

Klj  =  -o^S‘C(^i  +  v'^c  ^ 

^  ot  K  C '  V)  )  C 

=  o(.kC\+v)(a ^  c &  +  0 ) 

at  the  top  and  bottom  edges,  and  having  edge  deformations  given  by, 

~>t 

AaX 


—  —  A  ^s  ^ 


a 


it  is  clear  that  the  solution  presented  in  this  report  will  satisfy 
these  loading  and  edge  conditicais  and  will  conse<iaently  represent  the 
required  basic  conical  element  required  for  incorporation  in  an  extended 
version  of  SM-385OO. 
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